Accepted for Publication in the Astrophysical Journal 
Preprint typeset using LAT^^X style cmulateapj v. 20/04/00 



o 
o 

(N 

< 
O 



> 
O 
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ABSTRACT 

We investigate available constraints on the circumstellar medium (GSM) around long-duration 7- 
ray burst (GRB) progenitors from afterglow spectra. We first establish a statistical sample of five 
GRB afterglow spectra that have been collected and analyzed with no prior knowledge of the line-of- 
sight properties. This sample is then adopted for a uniform search of Wolf-Rayet wind signatures, as 
represented by CIV AA 1548, 1550 absorption doublets at Av = —1000 to —5000 km s~^ from the 
GRBs (hereafter CIV15). We report the detection of a single CIV15 absorber at Au « —1500 km s~^ 
from GRB 050730 and none in the rest. Our search yields an estimate of 20% for the incidence of 
CIV15 absorbers with rest-frame absorption equivalent width EW(CIV1548) > 0.2 A near GRB host 
galaxies, consistent with the incidence of intergalactic C IV15 near classical damped Lya absorbers toward 
quasar sightlines. Including the two CIV15 absorbers previously known toward GRB 021004, we further 
demonstrate that the presence of H°, C+, and Si+ together with the absence of excited C+ or Si"*" argue 
against a GSM origin. The null result is consistent with the expectation that the circumburst medium 
is fully ionized by the afterglow radiation field. We examine possible scenarios for the survival of the 
C'^+ ions, including a clumpy wind model. We find that a clumpy wind is unable to effectively shield the 
ionizing radiation and allow C'^"'" to survive at r < 10 pc from the afterglow. We conclude that the lack 
of CSM-originated C IV15 absorbers is consistent with Wolf-Rayet winds terminating at < 30 pc from 
their progenitor stars. 

Subject headings: gamma rays: bursts — ISM: abundances — ISM: kinematics — intergalactic medium 
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1. INTRODUCTION 

A massive star origin of most long-duration, soft-spectrum 
7-ray bursts (GRBs) has now been well-established on em- 
pirical grounds (see Wooslcy & Bloom 2006 for a recent 
review) . Not only the locations of these GRBs trace closely 
the spatial distribution of young stars in their host galax- 
ies (e.g. Bloom, Kulkarni, & Djorgovski 2002, Hammer et 
al. 2006; Fruchter et al. 2006), but a growing number of af- 
terglows are also observed to exhibit spectral features that 
are characteristic of a core-collapse supernova (Hjorth et 
al. 2003; Stanek et al. 2003; but see Fynbo et al. 2006a, 
Gal- Yam et al. 2006, & Delia Valle et al. 2006). In the 
popular coUapsar model, likely candidate progenitors of 
long-duration GRBs are Wolf-Rayet stars that undergo a 
significant mass loss through stellar winds (MacFadyen & 
Woosley 1999). It is therefore expected that the circum- 
stellar medium (GSM) has been regulated by the ejected 
wind material throughout the lifetime of the progenitor 
star (Wijers 2001; Ramirez-Ruiz et al. 2001; Chevalier et 
al. 2004; Ramirez-Ruiz et al. 2005; van Marie et al. 2005; 
Eldridge et al. 2006; van Marie et al. 2006). 

In principle, the presence of Wolf-Rayet winds in the 
circumburst medium can be demonstrated in an early- 
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time, single-epoch afterglow spectrum. Specifically, ul- 
traviolet spectra of local Wolf-Rayet stars exhibit promi- 
nent P-cygni profiles in N V A 1242, Si IV A 1442, and C IV 
A 1550 trannsitions (e.g. Niedzielski & Rochowicz 1994; 
Willis et al. 2004) that indicate a wind terminal velocity 
up to Uterm — 5000 km s"^ (Figure 1). Numerical simula- 
tions of single stars show that a free-expanding wind will 
imprint a narrow absorption-line component in the after- 
glow spectrum, blue-shifted (Aw < —1000 km s^^) from 
the host redshift (van Marie, Langer, & Garci'a-Segura 
2005). Such blue-shifted, high- velocity absorption compo- 
nents have been reported in C IV transitions for GRB 021004 
at 0GRB = 2.329 (M0ller et al. 2002; Mirabal et al. 2003; 
Schaefer et al. 2003; Fiore et al. 2005; Fynbo et al. 2005; 
Starling et al. 2005; see also Lazzati et al. 2006), GRB 030226 
(Klose et al. 2004), and GRB 050505 at zqrb = 4.2741 
(Berger et al. 2006), but the origin of these high- velocity 
C IV absorption components is ambiguous for three main 
reasons. 

First, over a large velocity interval, between Av — 1000 
and Av = 5000 km s~^ blueward of the host redshift, 
the likelihood of finding an intervening absorber origi- 
nated in a foreground galaxy along the line of sight is non- 
negligible. Second, the presence of an intense ultraviolet 
radiation field from the optical afterglow substantially in- 
creases the ionization of the circumburst medium. It there- 
fore prohibits most ions from surviving at < 30 pc radius 
from the afterglow. A GSM origin of these absorbers would 
require the Wolf-Rayet winds to reach beyond 30 pc, but 
observations of local Wolf-Rayet stars show that the wind 
blown bubbles have a typical size of < 10 pc (Gruendl et 
al. 2000). Finally, the circumburst medium is expected to 
be highly ionized even at > 30 pc, due to the presence 
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of the intense afterglow radiation field. But the observed 
CIV abundance alone is insufficient for constraining the 
ionization state of the gas around the afterglow. A full 
consideration of all species at different ionization stages is 
necessary (e.g. Fynbo et al. 2005). 

We have established a statistical sample of five GRB af- 
terglows for which the afterglow spectra were collected and 
analyzed with no prior knowledge of the absorption-line 
properties along the line of sight. We have also assembled 
an early sample of previous GRB afterglows for which the 
spectra are available to us. The statistical sample allows 
an unbiased estimate of the incidence of the blue-shifted, 
high- velocity CIV absorption components. The statistical 
sample and early sample together allow us to investigate 
the nature of the blue-shifted, high- velocity C IV absorp- 
tion components based on a suite of ionization models. In 
§ 2, we present the GRB samples and summarize the avail- 
able spectroscopic data for each burst. In § 3, we examine 
the possible intergalactic nature of these CIV absorbers 
based on the relative abundances between different ion- 
ization states of carbon and silicon. In § 4, we consider 
the afterglow radiation field estimated from light-curve ob- 
servations, and constrain the extent of Wolf-Rayet winds 
based on the presence and absence of high- velocity C IV 
absorbers in the afterglow spectra. We further discuss the 
scenarios that will allow C^+ to survive the intense after- 
glow radiation field. Finally, we summarize our study in § 
5. 

2. THE GRB SAMPLES 

We have compiled a statistical sample of GRB sight- 
lines for which we resolve the CIV AA 1548, 1550 absorp- 
tion doublet (i.e. a spectral resolution 6v < 180 km s~^) at 
signal-to- noise ratios oi S/N > 7 per resolution element so 
that the C IV A 1548 transition would be observed at > 3 cr 
confidence level to an EW limit of 0.4 A over a 170 km s^^ 
velocity interval. The GRB sightlines in this sample are 
selected entirely based on the data quality (spectral reso- 
lution S V and S/N) from the GRAASFQ database with no 
prior knowledge of the presence or absence of intervening 
absorbers. These include GRBs 050730, 050820, 050908, 
051111, and 060418, and form a statistically unbiased sam- 
ple for a blind search of intervening C IV absorbers. 

In addition, we examine the GRB spectra that are avail- 
able in the literature, in search of Wolf-Rayet signatures. 
These sightlines, despite incomplete and biased by known 
light-of-sight properties, offer additional insights for con- 
straining the Wolf-Rayet nature of GRB progenitor stars. 
The GRBs in this catagory are GRBs 000926, 021004, and 
030323. We call this an "early sample" , because the events 
occurred prior to 2004. 

In both the statistical and the early samples, we search 
for blue-shifted, high- velocity components of C IV within 
|At>| < Wtcrm km s^^ in the afterglow spectra. The veloc- 
ity threshold wterm is chosen to match the terminal velocity 
observed for galactic Wolf-Rayet stars (Figure 1). Based 
on 131 objects in the Vllth catalogue of galactic Wolf- 
Rayet stars by van der Hucht (2001), we show that all have 
a terminal wind velocity of < 5000 km s^^. The search of 
high- velocity C IV components is therefore carried out at 
\Av\ < Wtcrm — 5000 km s^^. An additional selection 
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threshold for Au = —1000 km s^^ is applied based on the 
maximum velocity of galactic winds observed in starburst 
galaxies at z > 1 (e.g. Adelberger et al. 2005). The nature 
of CIV within Av > —1000 km s~^ becomes more am- 
biguous due to possibility of a large-scale galactic wind. 
We therefore exclude this region from consideration. 

The blue-shifted, high- velocity CIV absorbers identi- 
fied over the velocity interval of —5000 < Av < —1000 
km s~^ (hereafter designated as CIV15) from the GRB 
host redshift 2;grb represent likely wind signatures of the 
progenitor stars. We also search for the presence of CII, 
Sill, Si IV, and additional ions whenever the transitions 
are covered. These additional transitions are crucial for 
constraining the ionizing radiation field around the ab- 
sorbing gas. Here we summarize available spectroscopic 
data and our search results. For each source, the zero 
relatively velocity Av = corresponds to the redshift of 
the afterglow, which is determined based on the absorp- 
tion profiles of low-ionization transitions such as C II, Si II, 
AlII, and Fell. 

2.1. The Statistical Sample 
2.1.1. GRB 050730 at z = 3.96855 

This GRB was first discovered by the Swift satellite 
on UT 2005 July 30 (Holland et al. 2005). An optical 
source was found promptly using the Ultraviolet-Optical 
Telescope (UVOT) on board Swift with V = 17.6 about 
three minutes after the burst trigger (Holland et al. 2005). 
This source was subsequently noted by Cobb & Bailyn 
(2005) to fade, confirming it as the optical afterglow. We 
obtained an echelle spectrum of the afterglow, using the 
MIKE spectrograph (Bernstein et al. 2003) on the Mag- 
ellan Clay telescope, four hours after the initial trigger. 
Descriptions of the data were presented in Chen et al. 
(2005) and Prochaska et al. (2006). The spectrum cov- 
ers a full spectral range from 3300 A through 9400 A with 
a spectral resolution oi Sv w 10 km s~^ at wavelength 
A = 4500 A and (5 w « 12 km s^^ at A = 8000 A. The 
host of the GRB exhibits a strong damped Lya absorp- 
tion feature with log N{}i I) = 22. 15 ±0.05, and abundant 
heavy ions with > 95% of the absorption confined within 
-200 < At; < 50 km s^^. 

Figure 2 shows the absorption-line profiles of C IV AA 1548, 1550, 
Si IV AA 1393, 1402, CH A 1334, and SiH A 1526 observed 
in the host ISM. The spectral region near the C IV absorp- 
tion doublet is contaminated by the atmosphere A-band 
absorption forest, but we identify a CIV15 absorber at 
Av = —1500 km s~^. Comparison with a quasar spec- 
trum obtained using the same instrument at the same 
site (shown as the magenta curve in Figure 2) confirms 
the lack of atmosphere absorption transition at the loca- 
tion of these CIV15 lines. The identification is further 
supported by a consistent kinematic profile of the Si IV 
AA 1393, 1402 doublet transitions, for which we measure 
log 7V(SiIV) = 13.9 ±0.1. For the CIV transitions, we 
measure a lower limit to the rest-frame absorption equiv- 
alent width of EW(A 1550) > 0.36 A. At this redshift, 
the CII A 1334 transition is contaminated by a pair of 
CIV absorption doublet at lower redshifts, z = 3.254 and 
z = 3.258. Therefore, no constraints can be obtained. We 
do, however, identify weak Si II A 1526 absorption and a 
strong Si II A 1260 profile. 
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2.1.2. GRB 050820A at z = 2.6U7 

This GRB was first discovered by the Swift satellite 
on UT 2005 August 20 (Page ct al. 2005). An optical 
transient, reported less than 1 hour after the GRB, was 
identified in data taken shortly after the trigger (Fox & 
Ccnko 2005; Vcstrand ct al. 2006). Wc obtained an echelle 
spectrum of the afterglow, using the HIRES (Vogt et al. 
1994) on the Keck I telescope, an hour after the initial trig- 
ger. Descriptions of the data were presented in Prochaska 
et al. (2006). The spectrum covers a full spectral range 
from 3800 A through 8000 A with a spectral resolution 
of (5 w « 7.5 km s^^ across the entire spectral range. The 
host of the GRB exhibits a strong damped Lya absorp- 
tion feature with log A'^(HI)=21.0±0.1, and abundant 
heavy ions with > 95% of the absorption confined within 
-250 < At; < 200 km s"^ 

Figure 3 shows the absorption-line profiles of C IV AA 1548, 
Si IV AA 1393, 1402, Oil A 1334, Si II A 1260, and Al II A 1670 
observed in the host ISM. While the host ISM shows com- 
plex, multiple-component absorption features at \Av\ < 
250 km s~^ in CIV and Si IV transitions, it does not ex- 
hibit additional blue-shifted, high-velocity component be- 
yond this velocity range at a 3-a upper limit of EW(A 1548)= 
0.2 A over 6v = 7.5 km 

2.1.3. GRB 050908 at z = 3.344 

This GRB was discovered by the Swift satellite on UT 
2005 September 8 (Goad et al. 2005). An optical tran- 
sient, seen in data taken 14 minutes after the burst, was 
reported « 80 minutes after the GRB trigger (Torii et al. 
2005; see also Cenko et al. 2005 and Li et al. 2005). We 
obtained a low-resolution spectrum of the afterglow, using 
CMOS on the Gemini north telescope, four hours after the 
initial trigger (Foley ct al. 2005). A moderate-resolution 
spectrum was taken four hours later by our group using 
DEIMOS (Faber et al. 2003) on the Keck II telescope 
(Prochaska et al. 2005). The final stacked CMOS spec- 
trum covers a spectral range from 5073 A through 7945 A 
wfth chip gaps at 6006-6022 A and 6972-6987 A a spec- 
tral resolution of i5v « 200 km s~^. The final stacked 
DEIMOS spectrum covers a spectral range from 6335 A 
through 8960 A with a spectral resolution of 5v w 40 
km s~^. 

The CMOS spectrum covers the redshifted hydrogen 
Lya A 1215 transition. We measure log 7V(II I) = 19.2 ± 
0.2. Figure 4 shows the absorption-line profiles of CIV 
AA 1548, 1550, Si IV AA 1393, 1402, C II A 1334, Si II A 1526, 
and AlII A 1670 observed in the host ISM. The host galaxy 
also exhibits abundant heavy ions with > 95% of the ab- 
sorption confined within —300 < Az; < 200 km s^^. We 
do not detect CIV15 beyond this velocity range at a 3-cr 
upper limit of EW(A 1548)= 0.28 A over ^ u = 40 km 3"^ 

2.1.4. GRB 051111 at z = 1.54948 

This GRB was discovered by the Swift satellite on UT 
2005 November 11 (Sakamoto et al. 2005). An optical 
transient was found 27 s after the burst on the ground 
(Rykoff et al. 2005). An echelle spectrum was taken an 
hour after the initial burst by the Keck Observatory staff, 
using HIRES on the Keck I telescope. However, the echelle 
spectrum does not extend below 4000 A, where the red- 
shifted CIV absorption is expected. We also obtained a 



low-resolution spectrum of the afterglow, using CMOS on 
the Gemini north telescope, 2.6 hours after the initial trig- 
ger. The final stacked GMOS spectrum covers a spectral 
range from 3800 A through 4400 A with a spectral resolu- 
tion of Sv ~ 170 km s""'^. Figure 5 shows the absorption 
profiles of C IV AA 1548, 1550, Si II A 1526, and Al II A 1670 
observed in the host ISM. The host galaxy also exhibits 
abundant heavy ions with > 95% of the absorption con- 
fined within -200 < A w < 100 km s'^ We do not detect 
CIV15 beyond this velocity range at a 3-(t upper limit of 
EW(A 1548) = 0.37 A oveidv = 170 km s'^. 

2.1.5. GRBO6O4I8 at z = 1.4901 

This GRB was first discovered by the Swift satellite on 
UT 2006 April 18 (Falcone et al. 2006). A candidate op- 
tical afterglow, with V = 14.5 mag, was reported from 
155UyOT images taken starting 88 seconds after the trigger 
(Falcone et al. 2006). We obtained an echelle spectrum 
of the afterglow, using the MIKE spectrograph on the 
Magellan Clay telescope, 28 minutes after the initial trig- 
ger. Descriptions of the data were presented in Prochaska 
et al. (2006). The spectrum covers a full spectral range 
from 3500 A through 9000 A with a spectral resolution of 
(5?; « 12 km s"^ at wavelength A = 7000 A. The host of 
the GRB exhibits abundant heavy ions with > 95% of the 
absorption confined within —200 < Av < 100 km s^^. 

Figure 6 shows absorption-line profiles of C IV AA 1548, 1550, 
Si IV AA 1393, 1402, C II A 1334, Si II A 1526 and Al II A 1670 
observed in the host ISM. Over the velocity interval from 
A f = —5000 km s^^ through Av = —1000 km s^^, we do 
not detect additional components at a 3-a upper limit of 
EW(A 1548) = 0.36 A over a spectral resolution element 
of (Ju = 12 km s~^. 

2.2. The Early Sam,ple 

We have searched CIV15 along the sightlines toward 
GRB 000926 (Hurley et al. 2000; Harrison et al. 2001), 
GRB 021004 (Shirasaki et al. 2002; Holland et al. 2003), 
and GRB 030323 (Graziani et al. 2003; Gilmore et al. 2003), 
for which afterglow spectra are available to us. The af- 
terglow spectrum of GRB 000926 was obtained using the 
Echelle Spectrograph and Imager (ESI; Shcinis et al. 2002) 
on the Keck II telescope (Castro et al. 2003). We found 
no CIV15 at a 3-cr upper limit to the rest-frame absorp- 
tion equivalent width EW(A 1548) = 0.35 A over a spec- 
tral resolution element of 6v = km s~^. The afterglow 
spectrum of GRB 030323 was obtained using the F0RS2 
spectrograph on the VLT/UT4 telescope (Vreeswijk et al. 
2004). We found no C IV15 at a 3-ct upper fimit to EW(A 1548) = 
0.04 A over a spectral resolution element of 5v = 143 
km s~^. 

Afterglow spectra of GRB 021004 were obtained, using 
various spectrographs on multiple telescopes (Castander 
et al. 2002; Chornock & Filippenko 2002; Salamanca et 
al. 2002; Castro-Tirado et al. 2002; MoUer et al. 2002; 
Fiore et al. 2005; Starling et al. 2005). The host of this 
GRB does not exhibit a prominent damped Lya absorp- 
tion feature that is commonly seen in most GRB hosts (e.g. 
Vreeswijk et al. 2004; Chen et al. 2005; Berger et al. 2006; 
Jakobsson et al. 2006). Based on the Lyman limit absorp- 
tion discontinuity and the observed absorption strength of 
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the Ly/3 transition, Fynbo et al. (2005) estimated a neu- 
tral hydrogen column density of log A^(HI) = 19.5±0.5. 
In addition, available echelle spectra of the optical after- 
glow clearly resolve two C IV15 absorption components ap- 
proaching Av — —3000 km s~^. While a blue-shifted, 
high-velocity component is expected in the circumstellar 
medium of a Wolf-Rayet progenitor, the presence of two 
components is difficult to explain under the free-expanding 
wind model of a single progenitor star (e.g. van Marie et 
al. 2005). 

We retrieved from the ESO Science Archive (program 
ID'S 7Q.A-0599(B) and 70.D-0523(D)) echelle spectra of 
the source obtained 13.5 hours after the burst, using the 
Ultraviolet- Visual Echelle Spectrograph (UVES; D'Odorico 
et al. 2000) on the VLT Kueyen telescope. We exam- 
ined and reduced the individual spectra. The data were 
weighted by their S /N ratios and co-added to form a final 
stacked spectrum. The final, stacked spectrum spanned a 
spectral range from 4200 to 9900 A with a spectral reso- 
lution of w 3.8 km s^^ per pixel, and achieved a signal-to- 
noise ratio of S/N « 20 per resolution element at around 
5200 A. A copy of a low-resolution spectrum that was ob- 
tained nearly 4 days after the GRB, using the Low Resolu- 
tion Imaging Spectrograph (LRIS) on the Keck I telescope 
was kindly made available to us by N. Mirabal. The spec- 
trum covers the spectral range around the redshifted Lya 
transition at the host redshift and has a spectral resolu- 
tion of w 150 km s^^ per pixel. Descriptions of the LRIS 
observations are presented in Mirabal et al. (2003). 

We present in Figure 7 the absorption-line profiles of 
HI Lya A 1215, CIV AA 1548, 1550, Si IV AA 1393, 1402, 
CII A 1334, Sill A 1526, and AlII A 1670 observed in the 
host of GRB 021004. Figure 7 shows that in addition to 
the main absorption component at |At;| < 400 km s^^ 
two absorption components are seen blue-shifted from the 
host redshift at |At;| = 2700 (Component 1) and 2900 
km (Component 2), respectively. A strong Lya ab- 
sorption feature is present at the location of the two blue- 
shifted, high-veolicty C IV components, but not resolved in 
the LRIS spectrum. Mirabal et al. reported a rest-frame 
absorption equivalent width of EW(Lya)= 3.9 ± 0.6 A^^ 
which implies a total neutral hydrogen column densitjcl 
of log iV(II I) ~ 19.5 for a typical Dopplcr parameter of 
6 = 20 km s~^. Our A^(H I) estimate is consistent with 
that reported by Fynbo et al. (2005) from their HST ob- 
servations of the Lyman limit discontinuity. We confirm 
this value with a Voigt profile fit to the LRIS data and 
estimate an uncertainty of 0.4 dex which is dominated by 
systematics (continuum fitting and line-blending). In ad- 
dition to the strong HI Lya transition, we also observe 
abundant low-ions such as C"*" as evident by the C II A 1334 
transition at the location of Component 1. 

We measure the ionic column density of each species us- 
ing the apparent optical depth method (Savage & Sembach 
1991). The measurements are presented in Table 1. The 

^Mirabal et al. (2003) derived progressively larger A'(H I), from 
log iV(H I) = 14.95 to log 7V(H I) = 16.11, based on the observed 
EW(Lya), EW(Ly/3), and EW(Ly7). These authors assumed that 
these transitions fall on the linear part of the curve of growth, but 
noted the differences in the derived A'^(H I) as suggestive of the lines 
being saturated. The differnce between our derived A'(H I) and those 
from Mirabal et al. arises entirely due to the assumption these au- 
thors adopted, which is invalid for saturated lines. 



large contrast in the observed relative column densities of 
C"*" and C'^^ for components 1 and 2 suggests that they 
arise in distinct phases (see § 3.2). Furthermore, strong ab- 
sorption from neutral hydrogen and low ionization species 
such as C"*" strongly indicate that this absorber may not 
originate in the highly ionized circumstellar medium of the 
progenitor. We will discuss various ionization scenarios in 
§3.2. 

3. THE NATURE OF BLUE-SHIFTED, HIGH- VELOCITY CIV 
ABSORBERS 

Our systematic search for C IV15 absorbers along five 
GRB lines of sight in the statistical sample has yielded 
identifications of one absorber along the sightline toward 
GRB 050730 at Aw « -1500 km s'^ Despite a het- 
erogeneous spectral data set, our survey is sensitive to 
C IV15 absorbers of rest-frame absorption equivalent width 
EW(A 1548) > 0.04 - 0.4 A at the 3-cr significant level, 
when scaled to a uniform Sv = 170 km s~^ velocity reso- 
lution element. We estimate that the incidence of CIV15 
absorbers with EW(A 1548) > 0.4 A is 20% near a GRB 
afterglow. 

Despite a small sample, we attempt to determine the 
confidence level of our estimate in the incidence of C IV15 
absorbers based on a bootstrap re-sampling analysis. We 
first adopt the five GRB sightlines as the parent sample. 
Then we randomly select five fields from the parent sam- 
ple, allowing for duplication. Next, we measure the inci- 
dence of C IV15 absorbers in the new sample. Finally, we 
repeat the procedure 1000 times and measure the scatter 
in our measurements. Our bootstrap analysis yields a 68% 
confidence interval of — 40% in the incidence of C IV15 
absorbers with EW(A 1548) > 0.4 A. 

Including sources in the early sample from the litera- 
ture (GRBs 000926, 021004, and 030323), we confirm addi- 
tional two components along the sightline toward GRB 021004 
at A w « —2800 km s~^. We note that possible C IV15 fea- 
tures have also been reported for GRB 030226 at zgrb = 
1.986 by Klose et al. (2004) and discussed by Shih et al. 
(2006). The high spectral resolution data are not accessi- 
ble to us. The source is therefore not included in our analy- 
sis. A summary of the observed absorption line properties 
is presented in Table 2, where we list for each field the 
estimated isotropic energy release of 7-ray photons Eiso, 
the redshift of the GRB host, A^(H I) and metallicity of 
the host ISM, the maximum velocity span of the observed 
blue-shifted CIV absorption doublet, and the time period 
when the afterglow spectra were obtained fj^^^. We find 
no correlation between the presence of C IV15 and other 
factors, such as Eiso or t^J^^'^. 

While a Wolf-Rayet wind from the GRB progenitor is a 
natural candidate for explaining the observed C IV15 ab- 
sorbers, it is clear that the majority of GRB afterglow 
spectra do not exhibit these features. Here we investigate 
the nature of these C IV15 clouds, taking into account the 
presence of ions at lower ionization stages. 

3.1. The Likelihood of Contaminations by Intervening 
Absorbers at Zabs < ^grb 

We first consider the statistical sample of five GRB lines 
of sight as a whole and evaluate the likelihood of finding 
an intervening CIV15 absorber over |At>| = 1000 — 5000 
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km s""'^ at zciv < zgrb based on known statistics from 
quasar absorption-line studies. A non-negligible proba- 
bility of finding a foreground absorber over the redshift 
pathlength surveyed by this sample would challenge the 
validity of attributing these absorbers to the GRB host 
environment. 

Random surveys toward quasar lines of sight have shown 
a mean number density per unit redshift interval of ncw{z = 
2) = 1.8 ± 0.4 for CIV absorbers of EW(A1548) > 0.4 A 
and ncw{z = 2) - 5 for EW(A1548) > 0.15 A (Steidel 
1992). The total redshift pathlength over Av = 1000 - 
5000 km s~^ along the GRB lines of sight in our sample of 
five is A z = 0.23. Therefore, the probability of detecting 
at least one random C IV absorber of EW(A 1548) > 0.4 
is 34%. At lower threshold EW(A1548) > 0.15, the prob- 
ability of detecting at least one random absorber at lower 
redshifts is 64%. We note that the C IV15 absorber toward 
GRB 050730 has EW(A 1548) > EW(A 1550) > 0.36 A. 

It is also necessary, however, to take into account the 
possibility of large-scale clustering between C IV absorbers 
through clustering of the absorbing galaxies. Galaxy clus- 
tering may be important in this analysis, particularly be- 
cause the search is carried out near a GRB host galaxy. 
Unfortunately, the cross-correlation function of interven- 
ing C IV absorbers and GRB host galaxies is not well de- 
termined. Existing two-point correlation measurements 
are for C IV absorbers of a wide range of observed column 
density (e.g. Petitjean & Bergeron 1994). These arc likely 
underestimates of the clustering extent and amplitude for 
strong CIV absorbers, because stronger absorbers are ex- 
pected to cluster more strongly with one another (Adel- 
berger et al. 2005). 

As a first estimate, we have searched a sample of 53 
quasar sightlines from the Keck/UCSD High Resolution 
database (Prochaska et al. 2007) with known intervening 
super Lyman limit systems (SLLS) of log iV(II I) > 19 
or damped Lya absorbers (DLAs) of log A''(H I) > 20.3 
(comparable to the N{li I) found in GRB host galaxies) 
at 1.7 < Zabs < 4.3. The DLAs are believed to trace the 
majority of the galactic population at z > 2. Although 
their luminosity (or mass) function has not been estab- 
lished empirically, current expectation is that a significant 
fraction correspond to sub-i* galaxies (e.g. Weatherley et 
al. 2005) similar to GRB host galaxies (e.g. Christensen, 
Hjorth, & Gorosabel 2004). 

We searched for C IV15 absorbers in the vicinity (—5000 kms 
Av < —1000 kms~^) of intervening DLAs toward back- 
ground quasars (Figure 8). This QSO-DLA sample serves 
as a control sample for estimating the statistical signifi- 
cance of possible "over-abundance" of C IV15 absorbers in 
the vicinity of GRB host galaxies. In the sample of 53 
DLAs identified along quasar sightlines, we find a total of 
six CIV15 absorbers to an EW limit of 0.2 A near five 
DLAs. Therefore, the incidence of CIV15 from a DLA is 
« 0.11. Namely, there is on average one CIV15 absorber 
in the vicinity of every nine DLAs. Adopting this expecta- 
tion, we estimate that the probability of detecting at least 
one C IV15 absorber near a GRB host galaxy in our sta- 
tistical sample of five GRBs is 42%. While the sample of 
GRB host galaxies is still small, our exercise shows that the 
probability for the CIV15 absorber toward GRB 050730 
to originate in a foreground galaxy at ^; < 2:grb is non- 



negligible. 

3.2. Constraints from, the Absence of Excited, Ions 

A nearly generic feature in the afterglow spectra of the 
ISM surrounding the GRB is the presence of the fine- 
structure absorption of C"*", Si"*", 0° and Fe"*" (e.g. Chen 
et al. 2005; Berger et al. 2006; Fynbo et al. 2006). While 
these lines were initially interpreted as signatures of a high 
density, circumstellar medium, Prochaska, Chen, & Bloom 
(2006; hereafter PCB06) demonstrated that indirect UV 
pumping would dominate the excitation if the gas is at 
r < 100 pc of the afterglow. Therefore, the presence or 
absence of absorption due to excited ions constrains the 
distance of the gas from the GRB event (see also Shin et 
al. 2006; Vreeswijk et al. 2006). 

We examine the velocity profiles of various ions associ- 
ated with the C IV15 absorbers from GRB 021004 and from 
GRB 050730 (Figure 9). We include the low-ionization 
transitions of CII A 1334 and SiIIA1260 for GRB 021004 
and GRB 050730, respectively, together with their associ- 
ated fine-structure transitions. Both C+ and Si+ have a 
J = 1/2 ground-state and a corresponding J = 3/2 ex- 
cited fine-structure level. It is clear that the ground-state 
is signficantly populated for Component 1 of GRB 021004 
and for the majority of the line-profile for GRB 050730. 
We set upper limits to the ratio of the excited to ground 
level / = A^j=3/2/Af7=i/2 of / < 0.05 for Component 1 
of GRB 021004 and / < 0.1 for GRB 050730. Although 
optical depth effects will generally imply / < 1 (PCB06), 
the very low / values require the gas lay at large distance 
from the GRB afterglow. 

A lower limit to the distance of the gas from the af- 
terglow can be obtained by calculating the total number 
of excitations from the afterglow prior to the onset of the 
spectroscopic observations, f^^'^. For both of these GRBs, 
we can neglect de-excitations from spontaneous decay be- 
cause the lifetimes of the excited levels are considerably 
longer than i^^g'^. Similarly, if the electron density of the 
gas were sufficient to coUisionally de-excitc the states, one 
would predict / > 0.1 (Silva & Viegas 2000). Adopting 
GRB 021004 as an example, we find from Figure 12 that 
the optical afterglow light curve is very complicated for 
this GRB (e.g. Pandey et al. 2003; Bersier et al. 2003). 
We therefore derive a spline-fit model of the empirical 
^iJc-band measurements reported by Holland et al. (2003) 
a^ Uemura et al. (2003). In addition, we assume that 
oc with /3 = — 1 noting that our calculation is not 
sensitive to this choice. We must also model the absorp- 
tion line-profile of the gas 'cloud' to determine the num- 
ber of photons absorbed prior to the onset of our spec- 
troscopic observations. A single absorber with column 
density A^(CII) = lO-'^^'^cm^^ and a Dopplcr parameter 
6 = 10 km s~^ is a conservative approximation to the ob- 
served line-profile (the equivalent width actually under- 
estimates the observed value). In the following, we only 
consider the C II AA1036 and 1334 transitions. 

The number of photons absorbed from ions in the ground- 
state of C"*" prior to spectroscopic observations is estimated 
by integrating over the light curve from ti = 1000 s to t^^^'^. 
We choose ti = 1000 s because the light curve is not well 
defined prior to this time and to give a more conservative 
estimate of the photon number. Integrating to = s 
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would likely increase the photon number by < 30%. We 
find a photon surface density at distance r from the after- 
glow, 

iVphot = 10i«(^^)' cm-2 . (1) 

While only a fraction of the excitations would be followed 
by a spontaneous decay to the J = 3/2 level in the ground 
term, a conservative estimate is 50%, because the upper 
level has a higher J value and correspondingly higher sta- 
tistical weight. To set a lower limit to the distance of the 
gas from the afterglow, we calculate the distance at which 
the photon surface density matches twice the column den- 
sity limit for the J = 3/2 level of C+ as observed from 
the UVES observations, iV(CII*) < W^^cir-'^. This im- 
plies the gas occurs at a distance rmin > 1 kpc from the 
GRB afterglow. We derive a similar value for the Si"*" gas 
associated with GRB 050730. 

The above argument is based primarily on the direct 
observations of the afterglow flux. No extrapolation is re- 
quired for the prompt emission or to higher energies than 
observed for the afterglow radiation field. Furthermore, 
the conclusions are based on empirical observations of the 
abundances and are largely independent of shielding sce- 
narios (which we will consider to explain the survival of 
C^+; § 4.3). This is understood by the fact that optical 
depth effects would yield excitation of C+ in a thin layer 
of the gas facing the afterglow, producing a non-negligible 
amount of C II* transitions. For example, in a layer of 
0+ ions with A^(C+) = lO^^cm"^ where the C II transi- 
tions are optically thin, excited ions are expected to 
dominate but are absent in our data. 

However, it is possible to design a scenario in which 
layers of high-column density gas absorb a large fraction 
of iVphot and are consequently photo-ionized by the after- 
glow just prior to to^g'^. This would require significant fine- 
tuning. Altogether, the absence of fine-structure absorp- 
tion lends strong support for the conclusion that the C and 
Si ions identified atvfv 2800 km s~^ and v w 1500 km s~^ 
toward GRB 021004 and GRB 050730, respectively, are not 
circumstellar material of the GRB progenitor. 

3.3. Are the Absorbers Photo-ionized by the Metagalactic 
Radiation Field? 

If the C IVi5 absorbers found along the sightline toward 
GRB 021004 originate in foreground clouds at zc iv < zq^b , 
then these gaseous clouds are expected to be photo-ionized 
predominantly by the ultraviolet background radiation or 
possibly a local starburst (Simcoe 2006). Observations of 
abundance ratios between different ionization stages of a 
given species can be applied to constrain the ionization 
state of the gas. 

We use the Cloudy software (Ferland et al. 1998; ver- 
sion 06.02) to calculate the expected population ratios 
between different ionization stages of carbon and silicon 
for clouds of plane parallel geometry and under photo- 
ionization equilibrium. We have assumed a metallicity of 
0.1 solar and log A^(H I) = 15, appropriate for interven- 
ing, strong C IV absorbers along quasar lines of sight (e.g. 
Cowie et al. 1995). The results are not sensitive to the 
adopted metallicity and N(H I). We also note that the as- 
sumed solar relative abundances are not relevant because 



we compare pairs of ions for the same elemental. Figure 
10 shows the abundance ratios of C"*" to C"^+ and Si"^ to 
Si'^"'" versus the ionization parameter U = 4>^/cn]i, which 
represents the number of ionizing photons available per 
hydrogen particle. 

Column density measurements presented in Table 1 for 
both high- velocity C IV components allow us to place lim- 
its on the ionization parameter. For component 1, CII 
A 1334 is saturated. We can therefore place an upper 
limit at log U = —2.9. For component 2, both CIV and 
Si IV are saturated. The (CII/CIV) and (Sill/SilV) ra- 
tios yield a consistent lower limit of log U > —2.1. We 
note that log U < —3 is relatively low, even for Lyman 
limit systems (e.g. Bergeron & Stasfnska 1986; Lopez et 
al. 1999; Chen & Prochaska 2000). On the other hand, 
the CII/CIV ratio observed for the component 2 is typ- 
ical of what is obseved in the intcrgalactic medium (e.g. 
Songaila 2006). 

The inferred limits for the ionization parameter U allow 
us to examine the gas properties of the absorbers. Adopt- 
ing a meta-galactic radiation field of J912 = 1.4 x 10~^^ 
erg s~^ cm~^ Hz"-'^ sr'^ measured at z 2 (Scott et al. 
2000), we derive a gas density of nn ^ 0.01 — 0.1 cm~'^ 
for log ?7 = — 2 to log U = —3. Comparisons between 
observations and expectations of photo-ionized gas in the 
intcrgalactic medium therefore demonstrate that the ob- 
served relative ionic abundances for the two C IV compo- 
nents can be naturally explained under a photo-ionization 
equilibrium with the meta-galactic radiation field. In par- 
ticular, the upper limit on the ionization parameter for 
component 1 indicates that a higher gas density and lower 
ionization intensity is favored. Our photo-ionization anal- 
ysis substantiates the fine-structure analysis from the pre- 
vious section in that the gas is located at large distance 
from the GRB afterglow. 

If we are to interpret the C IV15 absorption associated 
with GRB 021004 as intcrgalactic, then one might expect 
to identify a galaxy associated with this gas. Fynbo et al. 
(2005) find six faint sources in HST imaging in a 10"xl0" 
region around the host. This overdensity of faint objects is 
qualitatively consistent with the overabundance of Mgll 
absorbers at lower redshift along the sightline. Fiore et 
al. (2005) reported two groups of Mgll absorbers identi- 
fied at z = 1.380 and z = 1.602 along the sightline (see 
also Prochter et al. 2006). Figure 11 shows these optical 
images obtained using HST/ ACS with the F606W filter 
three days (left panel) and eight months (right panel) af- 
ter the trigger. The images are retrieved from the HST 
data archive (Program ID — 9405). While the afterglow 
was still bright in the first epoch image, we point out a 
faint companion at an angular distance AO = 0.3" to the 
host galaxy, which corresponds to 2.5 kpc projected dis- 
tance at z = 2.299. The host galaxy and the companion 
together within a 0.5"-radius aperture have a total bright- 
ness of AB(F606W) = 24.55 ± 0.05 (Fynbo et al. 2005). 
We estimate that the companion contributes roughly 1/10 
of the total fiux and therefore has AB(F606W) « 27.3. 

Recall that a Lya absorption line with EW = 3.9 ± 0.6 

''Our analysis differs from those presented in Fiore et al. (2005) 
in two important aspects. First, we adopt only lower limits for sat- 
urated transitions. Second, we compare the same elements, but in 
different ionization stages. This is necessary due to potential biases 
imposed by different chemical abundance patterns in the CSM. 



Chen et al. 



7 



A is also identified at the location of the two C IV15 com- 
ponents (§ 2.2), although unresolved. The presence of Lya 
absorber of log N(H I) = 19.5 presents a further challenge 
to the origin of the absorber arising in the host environ- 
ment. On the other hand, the Lya absorption strength is 
consistent with originating in the extended gaseous halo of 
a faint dwarf galaxy at lower redshift (Chen et al. 2001). 
No deep image of the field around GRB 050730 has been 
reported yet. 

4. IMPLICATIONS FOR THE CIRCUMSTELLAR MEDIUM OF 
GRB PROGENITOR STARS 

The survey results from a statistical sample show that 
the majority of afterglow spectra do not exhibit expected 
wind signatures from Wolf-Rayet progenitors, consistent 
with the expectation that the circumburst medium is highly 
ionized. In addition, the fine-structure and photo-ionization 
analysis show that the known C IV15 absorbers near GRB 
afterglows most likely originate in a foreground galaxy 
along the sightline. Here we take into account the known 
radiation field from afterglow light-curve observations. The 
presence of wind signatures in an afterglow spectrum de- 
termines the minimum distance of the absorbing ions to 
the afterglow. The absence of wind signatures determines 
the maximum distance reached by the wind prior to the 
GRB phase, if the absence is due to photo-ionization of 
the wind. Considering a sample of afterglow spectra to- 
gether allows us to estimate the typical extent of Wolf- 
Rayet winds. In the following discussion, we consider both 
smooth and clumpy media for allowing effective shielding 
of ionizing photons from the afterglow. 

4.1. The Extent of Wolf-Rayet Wind around the 
Progenitor Stars 

Figure 12 shows the optical light curves of four GRB af- 
terglows in our spectroscopic sample, for which early-time 
( from < 10 minutes after the initial burst) photometric 
measurements are available. All but one of these after- 
glows reached an initial brightness of i? < 15.5 mag and 
declined by 2 magnitudes in less than two hours (with the 
exception of GRB 021004). Even at i? = 17.5, the cor- 
responding intrinsic luminosity easily competes with the 
brightest quasars known at this redshift (see e.g. Fan et al. 
2001) . Despite a brief lifetime, the afterglow can easily ion- 
ize the circumstellar medium to large radii (see PCB06). 
The presence or absence of the CIV15 doublet features, 
coupled with the radiation field measured from afterglow 
observations, therefore, determines the extent of the Wolf- 
Rayet wind from the progenitor stars. 

To assess the impact of the afterglow over the circum- 
stellar medium, we estimate the number of ionizing pho- 
tons released to the surrounding medium prior to the onset 
of the spectroscopic observations. We first adopt the af- 
terglow of GRB 050820 as an example, given its well sam- 
pled optical light curve that shows both the prompt emis- 
sion at tobs < ^GRB + 5 min and afterglow radiation at 
later times (Vestrand et al. 2006). The transient nature 
of the afterglow indicates that accurate constraints on the 
ionization fraction of the surrounding medium require a 
time-dependent photo-ionization calculation (e.g. Perna & 
Lazzati 2002). 

For an approximation of the ionization fraction of C^^, 
however, we note that the recombination rate coefficient 



aR{C^+) - 10-"-5 cm3 sec-i (Nahar & Pradhan 1997) 
at a gas temperature of T ~ 10^~^ K and electron density 
rie ^ 10^ cm~^, typical of a compact HII region, yields a 
characteristic recombination time scale 

irec(C4+)^ lyear» tl^. (2) 

Since the timescale for observation (< a few hours in the 
frame of the medium external to the GRB outflow) is 
much smaller than the recombination time for each GRB 
event (Table 2), the ionization fraction of the circumstel- 
lar medium can be approximated by considering only the 
amount of ionizing photons available to ionize C^+. 

For the afterglow of GRB 050820, we estimate an isotropic 
release of total ionizing photon number, (f)y{hiy > 64 eV) ~ 
1.8 X lO^" over At ^ t^'' = 3240 s. Defining r^in as the 
radius of the ionization front where 99.99% of C'^"'" are ion- 
ized, we have n(C3+)/n(C"*+) « 10""^ in an optically thin 
regime or 



n(C4 
Therefore, 



n(C3+) , 
^ ' = IQ-^ = cxp 



-0^ap,,(C3+) 



47rr: 



</'7^p/.(C3+) 



1.1 X 1039 



1/2 



pc. 



(3) 



(4) 



For (Tph{C^+) = 6.6 X lO^^^ cm^, we find that 99.99% of 
C3+ within 33 pc would be ionized. The lack of C IV15 in 
the spectrum is therefore consistent with the winds termi- 
nating at < 30 pc and becoming fully ionized shortly after 
the GRB event. 

The afterglow spectrum of GRB 050820 presented in § 
2.1.2 is sensitive to CIV transitions of log 7V2(t(CIV) > 
13.7 at > 2-a significance level for a nominal Doppler pa- 
rameter b = 15 km s~^. Therefore, it does not exclude the 
possibility that Wolf-Rayet winds with log N{C IV) < 13.7 
may be present at r > rmin- The lack of wind signatures 
in the host of GRB 050820 implies that if such winds exist 
and contain a large amount of 0^+ (log iV(CIV) > 13.7) 
then it did not reach beyond 30 pc. 

The presence of the two C IV15 components along the 
sightline toward GRB 021004 have been considered in the 
literature as the most likely candidate for Wolf-Rayet winds 
in the circumburst environment of GRBs, although our 
analyses in § 3 show that at least one (Component 1 at 
Au = —2675 km s~^) is best-explained by a foreground 
cloud located at > 1 kpc from the afterglow and photo- 
ionized by the ultraviolet background radiation field. 

Allowing the possibility for Component 2 at A w = —2900 
km s^^ to arise in the CSM of the progenitor star, we es- 
timate the minimum radius for the C^^ ions to survive 
the afterglow radiation field. Adopting the total num- 
ber of ionizing photons estimated by Lazzati et al. (2006) 
from both the prompt emission and the afterglow phase 
(j)-y = 4 X 10^", we derive a minimum radius rmin ^ 50 pc 
according to Equation (4) and find that more than 99% 
of the C3+ ions are ionized within this radiutQ. Combin- 
ing the upper limit derived for the extent of Wolf-Rayet 

*This is a factor of three smaller than the distance of the C IV 
clouds derived by Lazzati et al. The difference is primarily due to a 
different ionization fraction adopted by these authors. Comparison 
between Equation (3) in the present paper and Equation (11) in 
Lazzati et al. shows that these authors have allowed the gas to be 



60% ionized when defining rmin- This leads to rj, 



150 pc. 
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wind around the progenitor star of GRB 050820 (< 30 pc) 
with rinin derived for GRB 021004, we find a maximum ex- 
tent of 30 — 50 pc for the Wolf-Rayct winds around GRB 
progenitor stars. 

4.2. Can the C^^ Ions Survive the Afterglow Radiation 
Field atr <30 pc? 

The large distance implied for the Wolf-Rayet winds of 
the progenitor of GRB 021004 from observations of Com- 
ponent 2 at Au = —2900 km s^"'^ is much larger than 
what is commonly observed for local Wolf-Rayet stars. 
The typical wind terminal shocks are found at r < 10 
pc (Garci'a-Segura et al. 1996; Gruendl et al. 2000). In 
an effort to understand the physical conditions that de- 
termine the location of the wind terminal shock around 
a GRB progenitor star, van Marie et al. (2006) explored 
the parameter space spanned by the ISM density of the 
GRB host, the wind velocity, and the mass- loss rate dur- 
ing the Wolf-Rayet phase of the progenitor. According to 
their calculation, for the wind to reach beyond 50 pc at a 
wind-speed of t'wR = 3000 km only a narrow param- 
eter space is allowed. Their Figure 12 (van Marie et al. 
2006) shows that a mass-loss rate of 1 x 10^^ M© yr"-'^, an 
ISM density log pisM < —24 gcm~'^ or < 1 cm~'^ need to 
be satisfied. A lower ISM density or higher mass-loss rate 
allows the wind to reach large distances. 

A simple analytical treatment can give some useful in- 
sight. This is because a strict upper limit on the radial ex- 
tent of the Wolf-Rayet wind can be obtained by computing 
the radius of the termination shock for a free expanding 
wind. The radius at the inner edge of the wind bubble 
can be found by balancing the wind ram pressure with 
the post-shock cavity pressure. For a star that loses mass 
at a rate 10~^ M-s Mq yr~^ with a wind velocity 10^ Vw,3 
km in interstellar gas with density 10° nisM,o cm~^, we 
have an inner termination shock radius 



n = 15 



.V3/10 1/10 -3/10 2/5 



pc. 



(5) 



where 10^ tg is the lifetime of the star in Myr. For sufh- 
ciently high M or low nisM, the termination shock radius 
given in equation Q can become larger than Tmin- This 
requires 



A>3/10 1/10 -3/10 > o 
■'"-5 "w,3 "ISM.O ^ 



50 pc 



(6) 



which shows that if the wind is specially massive (M > 
10^* Mq yr^^) or the surrounding pressure is low (nisM < 
10^^ cm^"^), rt falls above the radial range for which 99.99% 
of the C'^"'" ions are ionized. Interestingly, this low gas den- 
sity is consistent with the low iV(II I) observed in the host 
galaxy of GRB 021004, log iV(H I) = 19.5 ± 0.5. 

Here we consider possible scenarios that would help to 
alleviate the burden of forcing the wind out to beyond 30 
pc. 

We first consider the structured jet model discussed in 
Starling et al. (2005). These authors proposed that the 
ionization of circumstellar medium is dominated by a nar- 
row and fast-moving jet, but the line of sight encompasses 

We note that at this ionization level the spectra cannot rule out the 
possibility that the observed A'^(CIV) is due to gas at r < 150 pc, as 
it would only imply a total gas column for carbon of log N{C) ~ 15 
within this radius. It is therefore necessary to adopt a more stringent 
ionization threshold. 



a wider and slower-moving jet through the wind, where 
the C'^+ ions can survive at smaller distances to the burst. 
We consider this scenario unlikely for two main reasons. 
First, following this structured jet model, we expect to ob- 
serve partial coverage of the absorption transitions along 
the line of sight. The CIV A 1548 absorption line profile 
of Component 2 is completely saturated over the central 
A Aobs — 0.8 A (Figure 7). We constrain the covering fac- 
tor of C'^^ ions along the line of sight at / > 97% at the 
2-a significance level. Second, light curve observations and 
modelling of GRBs 021004 and 050820 show that prompt 
emission of the GRB event contributes only a small frac- 
tion to the total ionizing photons (Vestrand et al. 2006; 
Lazzati et al. 2006). Therefore, the afterglow serves both 
as the background light probe and as the principle source 
of ionization of the circumstellar medium. 

Next, we consider the possibility of a clumpy wind for a 
chemical composition that is characteristic of local Wolf- 
Rayet winds. Wolf-Rayet type stellar winds are found to 
have significant inhomogeneous spatial distributions and 
a chemical composition dominated by helium and carbon 
with little hydrogen (Crowther 2007). Dense clumps, lo- 
cated at smaller distances to the afterglow along the line of 
sight can in principle serve as a natural screen for shield- 
ing the clumps that are further away from the progenitor 
stars. The same model is commonly applied to study the 
presence of dusty tori in the vicinity of quasars (e.g. Konigl 
& Kartje 1994). 

To test whether a clumpy wind is a viable solution, we 
use the compressible turbulent hot-star outflow descrip- 
tion of Moffat et al. (1994) as a working model, in which 
a simple scaling relation between the size, Ic, and den- 
sity, Pc of the clumps is assumed to prevail: pc oc l^^. 
What is more, Moffat et al. (1994) found that for Wolf- 
Rayet winds, a typical mass spectrum that follows to^^ ''^ 
provides a reasonably good description of the emission line 
variations, where mc oc pc Ic is the mass of a clump. Taken 
together, the above scaling laws favor an asperous wind 
with clumps of varying size but constant column density 
and a mass spectrum dramatically raising to lower clump 
masses. 

If full-scale turbulence really does prevail, the discrete 
clumps could in principle be usable to suppress a substan- 
tial fraction of the ionizing photon flux at Ej = hv > 64, 
eV from the afterglow. To calculate the properties of 
the asperous wind, we flrst consider the following three 
constraints. First, the clump column density should be 
at least 7V(CIV) 1.5 x lO^^ cm'^ so that 7V(CIV) x 
o'p/i(C'^+) = 1 for suppressing a substantial fraction of the 
ionizing ffux. Because the ionization potentials of He"*" and 
C'^^ are very similar (see the inset of Figure 13)0, we treat 
IIe+ and C'^"'" as the same particles in the photo-ionization 
process. This leads to a minimum gas clump gas column 
density of Ng = 7V(CIV) x (C/He)-i « 1.5 x lO^^ cm'^ 
for Wolf-Rayet winds with C/He — 0.1. Second, the clump 

°Wc note that photons of Epi^^t > 350 eV are energetic enough to 
remove electrons from the inner shell (K-shell states) and electrons 
in the outer shell (L-shell states) cascade down to fill the vacancy 
on time scales of a few fs (Schlachter et al. 2004). This explains the 
edge that appears at 350 eV in the inset of Figure 13. But because 
the afterglow radiation declines quickly with frequency (/3 = —1), 
ionization of due to the removal of a K-shell electron is limited 
to ~ 10%. 
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size should be at least comparable to the apparent image 
size of the afterglow (e.g. Granot et al. 2005), 

= 2 X IQie Eli% v^f [*days/(l + z)f" cm. (7) 
This leads to a miniiniim clump mean gas density of Pc.min ^ 
750 cm""^ and a minimum clump mass of mc,min 4 x 
1O~^M0. Third, the total mass of the unshocked Wolf- 
Rayet wind sets a limit to the total number of available 
clumps at the smallest scale (note that given the assumed 
mass spectrum, the total mass of the outflow is dominated 
by the lower mass clumps, here assumed to be ^ r±). In 
the unshocked wind, the mass within a radius r is Mr /vw, 
which combined with the minimum clump mass estimate 
gives the maximum number of clumps of size Ic ~ r± 
within the l/r^ wind: 

n. = W'M%'v-TPc,Lnr,c. (8) 
The volume filling factor is thus given by 

= 3 X 10-5 Af_5 vZl, p-^in ■ (9) 
Whether a clumpy wind scenario is a viable solution for 
the survival of at < 30 pc depends on the number of 
clumps one can reasonably find along the line of sight to 
the afterglow. Adopting a total ionizing photon emission 
of = 4 X 10^" for GRB 021004 from Lazzati et al. (2006), 
we estimate a photon flux accumulated in time of A^-y = 
3 X 10^^ cm~^ at r = 1 pc. For Pcmin one requires roughly 

n, = 2xl03p-i,i„r-2 (10) 

clumps to be present along the sightline to effectively ab- 
sorb all the ionizing photons. It is therefore clear that 
a clumpy wind model fails, because there is simply not 
enough mass in the unshocked wind medium. A close in- 
spection of the above requirements shows that for a fixed 
unshocked wind mass, altering Pc,min) T-pc or Ic (here as- 
sumed to be ~ r±_) fails to improve the situation. For 
example, in an extreme case scenario a single clump of 
Ng = 1.5 X 10^2 cm"2 (total clump mass of 0.05 M© 
for Ic = r±), at r=l pc along the sightline is sufficient 
to absorb the large incident ionizing flux, but requires 
M ^ lO~^M0yr~^ so that the probability of having one 
in front of the afterglow will be at least modest. We con- 
tend that a clumpy wind model is therefore an unlikely 
solution. 

After a careful consideration of various scenario, we con- 
clude that 0^+ ions at r < 30 pc cannot survive the intense 
afterglow radiation field to impose an absorption feature 
in the afterglow spectrum. 

5. SUMMARY AND CONCLUSIONS 

We have collected a sample of five GRB afterglow spec- 
tra for an unbiased search of blue-shifted, high-velocity 
CIV absorbers at |Ai;| = 1000 - 5000 km s'^ from the 
host redshift .zgrb- These absorbers (designated CIV15) 
are candidate stellar wind features expected from a Wolf- 
Rayet progenitor of long-duration GRBs. Our sample, al- 
though small, is statistically unbiased, because it is as- 
sembled based on all available afterglow spectra suitable 
for our study without prior knowledge of the line-of-sight 
properties. We identify only one C IV15 absorber at A t; = 
— 1500 km s^^ from GRB 050730 and none along the rest 
of the sightlines to a 3-<j limit of EW= 0.4 A over a uni- 
form 5v = 170 km s~^ velocity resolution element. 



Our search yields an estimate of 20% for the incidence 
of CIV15 absorbers from the GRB host galaxies with a 
68% confidence interval of 0-40%. This is consistent with 
what is expected for classical damped Lya absorbers to- 
ward quasar sightlines. The result suggests that the ma- 
jority of these C IV15 absorbers originate in a foreground 
galaxy along the sightline. 

We have also assembled an early sample based on GRB 
sightlines from earlier work, including GRBs 000926 (Cas- 
tro et al. 2003), 021004 (Mirabal et al. 2003; Fiore et al. 
2005), and 030323 (Vreeswijk et al. 2004). The sightline 
toward GRB 021004 is the only source that exhibits two 
CIV15 absorbers at Au = —2675 km (Component 1) 
and Aw = —2900 km s^^ (Component 2), respectively, 
as reported by previous authors. The CIV15 absorbers 
found in the early sample offer additional insights for un- 
derstanding their origin. 

The presence of a saturated C II A 1334 absorption and 
absence of C II* A 1335 transition together allow us to 
rule out the possibility of Component 1 originating in the 
GRB progenitor environment. Our analysis shows that 
had the gas originated in the vicinity of the afterglow (i.e. 
at r < 1 kpc) , it would have been highly excited by the 
UV photons of the afterglow. Instead, Component 1 most 
likely originates in a foreground galaxy along the line of 
sight, which is supported by both the statistical expecta- 
tion of the the incidence of C IV absorbers and the presence 
of a faint galaxy at 0.3" away from the GRB host galaxy. 

The lack of wind signatures for 80% of our statistical 
sample is understood as due to gas being photo-ionzed 
by the intense UV radiation field of the afterglows. Our 
estimates show that ionizing photons from the afterglow 
can ionize C'^+ to beyond 30 pc radius from the progeni- 
tor stars. We explore different scenarios that would allow 
C^+ to survive at close distances from the afterglow. We 
find based on a chemical composition C/Hc=0.1, typical 
of local Wolf-Rayet winds, a large density contrast in the 
wind can in principle provide necessary shielding for the 
0"^+ ions, but but requires M ^ 10~^MQyr~-'^ so that the 
probability of having one in front of the afterglow will be 
at least modest. We contend that a clumpy wind model is 
unlikely to be a viable scenario. 
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NNG05GF55G. H.-W.C. was partially supported by NASA 
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Table 1 

Summary of Absorption-line Properties 



Field 


Biso(xlO^=' ergs)» 


2GRB 


log N(H I) 


[fliSM 


" IV 




The Statistical Sample 


GRB 050730 ... 


1.2 ± 0.2 


3.968 


22.15 ± 0.05 


-2.0 ± 0.1 


1500" 


(4.0 - 5.7) 


GRB050820A . 
GRB 050908 ... 
GRB 051111 ... 


2.0 ± 0.4 
0.11 ± 0.01 
0.62 ± 0.05 


2.615 
3.344 
1.549 


21.00 ± 0.10 
19.2 ± 0.2 
> 20.3 


-0.6 ± 0.1 
~ -0.9 
< 0.9 


250 
200 
200 


(0.9 - 1.4) 
(3.7-4.3)7(7.7-8.4)'= 
(2.6 - 3.3) 


GRB 060418 ... 


0.67 ±0.05 


1.490 


> 20.3 


< 0.0 


250 


(0.5 - 2.3) 








The Early Sam 


iple 






GRB 000926 ... 


2.4 ±0.8 


2.038 


21.3 ± 0.2 


-0.17 ±0.05 


250 


(29.1 - 30.5) 


GRB 021004 ... 


0.5 ± 0.1 


2.328 


19.6 ± 0.3 


-1.3 ± 0.5 


2900 


(13.6 - 17.1) 


GRB 030323 ... 


0.3 ± 0.1 


3.372 


21.90 ±0.07 


-1.0 ± 0.2 


200 


(57.2 - 61.2) 



°^GRBs 050730, 050820A, 050908, 051111 & 060418: Butler et al. 2006, in preparation. GRB 000926: Harrison et al. 
2001; GRBs 021004 & 030323: Sakamoto et al. 2005. 

^Contaminated by the atmosphere A-band absorption. 

"^Spectra obtained using GMOS on the Gemini North telescope and DEIMOS on the Kecit telescope, respectively. See 
descriptions in § 2.6. 



Table 2 

Ionic C'oLrMx I)e\sitii;s of ■•llK:n-\-i:LOCTTY" Gas Toward GRB 021004''''' 



Transition 


component 1 
Av = -2675 km s'^ 


component 2 
At; = -2900 km s"^ 


CIV 1550 


14.07 ±0.02 


> 14.72 


CII 1334 


> 14.46 


13.62 ±0.04 


Si IV 1402 


contaminated 


13.8 ±0.3 


Sill 1526 


13.42 ± 0.04 


< 12.92 



'^We note that Fiorc et al. (2005) presented independent measure- 
ments for these transitions. However, they did not acc;oiuit for the 
saturation of strong lines and presented lower limits as in measure- 
ments with errors. 

'^Errors presented here account for only statistical errors. Uncer- 
tainties due to continuum fitting contribute to ~ 0.05 dex. 
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Fig. 1. — Number counts of Galactic Wolf-Rayet stars versus terminal wind velocity as determined from P-Cygni profiles of absorption lines 
like CIV A 1550 (van der Hucht 2001). 
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Fig. 2.— Absorption profiles of C IV AA1548, 1550, Si IV AA1393, 1402, CII A 1334, and Sill A 1526 observed in the host of GRB 050730. 
Zero relative velocity corresponds to z = 3.96855. Over the velocity interval from Av = —5000 km s~^ through Av = —1000 km s~^, we 
detect one blue-shifted CIV and SilV absorber at » -1500 km s'^ The CIV features are contaminated by the atmosphere A-band 
absorption forest. Although we are unable to obtain an accurate measurement of the column density of the ions, we place a lower 

limit to the rest-frame absorption equivalent width of EW(A 1550) > 0.36 A. The detection of the CIV absorption doublet is robust from 
the comparison with a quasar spectrum (the magenta curve) obtained using the same instrument at the same observation site (PKS2000; 
Prochter et al. 2007, in preparation). For SilV, we measure log A'^(SilV) = 13.9±0.1. The contaminating absorption features are dotted out. 
The dashed-dotted lines indicate the normalized continuum and zero level for guide. The 1-cr error array is shown in thin, cyan line. Zero 
relative velocity corresponds to z = 2.32897. 
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Fig. 3.— Absorption profiles of CIV AA1548, 1550, Si IV AA1393, 1402, CII A 1334, Sill A 1260, and AlII A 1670 observed in the host of 
GRB 050820. Zero relative velocity corresponds to z = 2.6147. Over the velocity interval from Av = —5000 km through Av = —1000 
km s~^, we do not detect additional CIV features with a 3-cr upper limit in the rest-frame absorption equivalent width of EW= 0.2 A over 
the spectral resolution element of Sv = 7.5 km s~^. 
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Fig. 4.— Absorption profiles of CIV AA 1548, 1550, Si IV AA 1393, 1402, CII A 1334, Sill A 1526, and AlII A 1670 observed in the host of 
GRB 050908. Zero relative velocity corresponds to z = 3.344. High-resolution data were taken using DEIMOS on the Keck II telescope and 
low-resolution data were taken using CMOS on the Gemini north telescope, which had non-contiguous spectral coverage due to chip gaps in 
the camera. One of the chip gaps is apparent in the Si IV doublet transitions. Over the velocity interval from Av = —5000 km s~^ through 
Av = —1000 km s~^, we do not detect additional CIV features with a 3-a upper limit in the rest-frame absorption equivalent width of 
EW(A 1548) = 0.28 A over the spectral resolution element of (Su = 40 km s~^. 
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Fig. 5.— Absorption profiles of CIV AA 1548, 1550, Sill A 1526, and AlII A 1670 observed in the host of GRB 051111 with CMOS on 
Gemini-N. Zero relative velocity corresponds to z = 1.54948. Over the velocity interval from Av = —5000 km s~^ through Av = —1000 
km , we do not detect additional CIV features with a 3-cr upper limit in the rest-frame absorption equivalent width of EW(A 1548) = 0.37 
A over the spectral resolution element of (5 = 170 km s"'^. 
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Fig. 6.— Absorption profiles of CIV AA 1548, 1550, Si IV AA 1393, 1402, CII A 1334, Sill A 1526, and AlII A 1670 observed in the host of 
GRB 060418. Zero relative velocity corresponds to z = 1.4901. Over the velocity interval from Av = —5000 km s~^ through Av = —1000 
km s~^, we do not detect additional CIV features with a 3-a upper limit in the rest-frame absorption equivalent width of EW(A 1548 = 0.36 
A over the spectral resolution element of (5 = 12 km s~^. 
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Fig. 7.— Absorption profiles of HI Lya A 1215, CIV AA 1548, 1550, Si IV AA 1393, 1402, CII A 1334, Sill A 1526, and AlII A 1670 observed 
in the host of GRB 021004. The same as Figure 2, where contaminating absorption features are dotted out. The dashed-dotted fines indicate 
the normaUzed continuum and zero level for guide. The l-cr error array is shown in thin, cyan line. Zero relative velocity corresponds to 
z = 2.32897. Over the velocity interval from An = —5000 km s~^ through Av = —1000km s~^, we identify two absorption components at 
Av = —2675km s~^ and Av = —2900km s~^, respectively. 
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Fig. 8. — A sample of four strong intervening hja absorbers, two DLAs and two super Lyman limit systems, at 2.2 < z < 3.2 toward 
quasar lines of sight. Contaminating features have been dotted out. These are part of the 53 quasar sightlines from the Keck/UCSD High 
Resolution database (Prochaska et al. 2007). We adopt this sample of 53 quasars as a control sample for estimating the statistical significance 
of possible "over- abundance" of C IV15 absorbers. 
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Fig. 9. — Velocity profiles for the CIV15 absorbers from GRB 021004 (left panels) and from GRB 050730 (right panels). We also present the 
low-ion transitions of CII A1334 and SiIIA1260 for GRB 021004 and GRB 050730, respectively, together with their associated fine-structure 
transitions. No excited transitions are detected in either of the sources. 
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Fig. 10. — Expected ionic abundance ratios versus the ionization parameter U = cfi-y/cnu. The expectations are calculated using the Cloudy 
software (Ferland et al. 1998; version 06.02) for clouds of piano parallel geometry and under photo-ionization equilibrium. We have assumed 
a metallicity of 0.1 solar and log A'^(H I) = 15 and adopted the background radiation field of Haardt & Madau (2006 in preparation) that 
includes ionizing photons from both quasars and starburst galaxies. Constraints on the ionization parameter based on the observed relative 
abundances between different ionization stages are shown in vertical lines for the two components toward GRB 021004. 
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Fig. 11.— Images of the field surrounding GRB 021004, obtained using HST/ACS with the F606W filter three days (left panel) and eight 
months (right panel) after the trigger. The second epoch image (right panel), obtained when the optical afterglow had faded, clearly shows a 
faint companion (pointed by the arrow) next to the host galaxies at the center of the frame. 



Chen et al. 



23 



."'■1 


1 


1 1 1 11^ — i—n 


W] 1 II III! 






N 




r 














iiiiJ 


J 


J III 


■1 1 



HP ID* 10^ JO^ 

nmj aiiMA HITB Trlj^ (^^^) 



10^ r 



I — I 1 1 111^ 



IQO LilJ I J H J I *im 

ID^ 10^ JO* Itf" HJ* 




Ifl" H)* 1^ 10* 10* 




IPP lOP 10^ W Ifl* 



Fig. 12.— Optical light curves of GRBs 021004, 050730, 050820, and 060418, for which photometric measurements from < 10 minutes after 
the initial burst are available. The vertical shaded regions indicate the period when the spectroscopic data were taJcen. The light curves of 

GRBs 021004 and 050730 arc based on R-band observations alone, while the light curves of GRBs 050820 and 060418 include observations 
iu the B, R, and I bands. We have corrected for the Galactic extinction, when combining photometry from different bandpasses with 
E{B -V) = 0.04 toward GRB 050820 and E{B -V) = 0.22 toward GRB 060418. Observations of GRB021004 were taken from Holland et al. 
(2003) and Uemura et al. (2003). Observations of GRB 050730 were taken from Klotz et al. (2005) and Damerdji et al. (2005). Observations 
of GRB 050820 were taken from Vestrand et al. (2006). Observations of GRB 060418 were taken from Schady & Falcone (2006) and Cobb 
(2006). 
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Fig. 13. — Estimated Af (C IV) versus distance to the afterglow in the presence of He+ for two different density profiles. The inset shows 
the photo- ionization cross sections of He+ and versus photon energy. The curves are calculated according to Verner & Yakovlev (1995) 
and Verner et al. (1996). 



